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The 7/ propagator in quenched QCD* 
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The calculation of the rj hairpin diagram is carried out in the modified quenched approximation (MQA) in 
which the lattice artifact which causes exceptional configurations is removed by shifting observed poles at k < k c 
in the quark propagators to the critical value of hopping parameter. By this method, the rj propagator can 
be accurately calculated even for very light quark mass. A determination of the topological susceptibility for 
quenched QCD is also obtained, using the fermionic method of Smit and Vink to calculate winding numbers. 
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1. INTRODUCTION 

The properties of the flavor singlet pseu- 
doscalar rj meson provide a unique phenomeno- 
logical window on the topological structure of 
QCD. The mass of the rf is believed to arise 
from topologically nontrivial gauge configurations 
via the axial U(1)a anomaly. In a semiclassi- 
cal treatment, the effect of instanton contribu- 
tions to the qq annihilation ("hairpin") diagram 
breaks the degeneracy between the rj and the fla- 
vor octet Goldstone bosons. In an effective chi- 
ral lagrangian description of QCD, the qq hair- 
pin can be interpreted as an 77' mass insertion. 
Within the large N c approximation, Witten and 
Veneziano [0 showed that the topological theory 
of U(1)a breaking is consistent with the mass- 
insertion view, and that in this approximation, 
the 7/ mass in the chiral limit is proportional 
to the topological susceptibility Xt of pure gauge 
(quenched) QCD. 
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(Note: Here we take f„ normalized to have the 
experimental value of m 96 MeV. This differs from 
that used in Ref. ||] by a factor y/2.) 

The study of the rj' propagator in lattice QCD 
is of great interest, not only as a quantitative 



check of the theory of the 17(1) problem, but also 
as a particularly sensitive probe of the differences 
between quenched and full QCD. For example, if 
the hairpin diagram behaves like a mass insertion, 
the quenched rj 1 momentum-space propagator is 
expected to include a double Goldstone pole con- 
tribution. 
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*Talk presented by H. Thacker 



As a result, quenched chiral logs arising from vir- 
tual 77' loops will complicate the chiral behavior 
of quenched QCD compared with that of the full 
theory ||. One of the purposes of the study re- 
ported here is to investigate in detail the time- 
dependence of the hairpin contribution to the 77' 
propagator and compare it with that expected 
from the double pole structure @. 

The method we use to calculate closed loops 
which originate at a given site of the lattice was 
introduced by Kuramashi, et al Q in their origi- 
nal study of the 7/ mass. In this "allsource" tech- 
nique, the quark propagator is calculated using 
a source given by a unit color-spin vector on ev- 
ery space-time point of the lattice. The closed 
quark loop from a given point is then calculated 
by assuming random-phase cancellation of other 
non-gauge-invariant terms. 

A particular difficulty encountered in the rj 
hairpin calculation is the presence of rapidly in- 
creasing non-gaussian errors in the limit of small 
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Figure 1. The unimproved hairpin propagator on 
a 12 3 x 24 lattice at = 5.7 and k = .1675 
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Figure 2. The MQA improved hairpin propagator 
on a 12 3 x 24 lattice at = 5.7 and k = .1675 



quark mass due to exceptional configurations. 
This problem is more severe in the hairpin cal- 
culation than it is in ordinary hadron spectrum 
calculations. Recently, the origin of the excep- 
tional configuration problem has been traced to 
the presence of topological zero mode poles in the 
quark propagator which have been displaced to 
values of the hopping parameter below k c by lat- 
tice effects. || A practical method for removing 
this lattice artifact, the modified quenched ap- 
proximation (MQA), has been proposed in Ref. 
Q . An example of the MQA improvement of the 
hairpin propagator is shown in Figs. 1 and 2. 
The propagators shown in Figs. 1 and 2 were ob- 
tained on a 12 3 x 24 lattice at = 5.7 with a 
hopping parameter of k = .1675 (m q ps 38 MeV). 
The improvement of errors due to the MQA pole- 
shifting is impressive. For lighter quark masses, 
the data for the 12 3 x 24 lattice without the MQA 
improvement is unusable due to extremely large 
errors. After MQA improvement it is possible 
to accurately calculate the hairpin diagram down 
to very light quark masses. In our calculations, 
we have shifted all poles which were found below 
k = .1690 and have calculated the hairpin for k 
values up to .1685. It may be feasible to go to 
even lighter quark mass values, but the calcula- 
tion of pole residues above k = .1690 becomes 
rapidly more expensive in computer time. 



2. THE rf PROPAGATOR 

The hairpin diagram has been calculated on 
two ensembles of gauge configurations available 
in the ACPMAPS library. One ensemble included 
200 configurations on a 12 3 x 24 lattice at = 5.7. 
The other ensemble consists of 200 configurations 
at = 5.7 on a 16 3 x 32 lattice. All together we 
have calculated the results for seven different val- 
ues of hopping parameter ranging from .161 to 
.1685. 

In general, the hairpin insertion may not be 
a simple p 2 -independent mass insertion, but in- 
stead may exhibit some p 2 -dependence. Expand- 
ing around the pion mass shell, we may write 

ml -> h(p 2 ) = ml + a(p 2 - jtvO + . . . (3) 

The second term corresponds to an additon to 
the 77' kinetic term in the chiral lagrangian. The 
form of the hairpin insertion may be determined 
by studying the time-dependence of the hairpin 
propagator. If the hairpin vertex is a simple p 2 - 
independent mass insertion, the measured propa- 
gator at zero 3-momentum should behave accord- 
ing to a pure double-pole formula, 

C 

A h (p = 0,t) = 2—3(1 + m^t)exp(-m w t) (4) 

+(*«-» (T-t)) 
The value of m w is determined quite accurately 
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Figure 3. The MQA improved hairpin propagator 
on a 16 3 x 32 lattice at j3 = 5.7 and k = .1680. 
The solid line is the time dependence of a pure 
double Goldstone pole, Eq. (4) expected if the rf 
hairpin vertex is a simple mass insertion. 



from the valence pion propagator, so the only ad- 
justable parameter in this fit is the overall nor- 
malization C. The second term in (3) would con- 
tribute a single-pole term to the propagator (2), 
which gives a term in (4) with pure exponential 
time dependence. Thus, both single and double 
pole terms are included by replacing the factor 
(7(1 + m n T) in (4) by (C 1 + C 2 m 7r t). The prelim- 
nary results of this analysis indicate that there is 
very little p 2 dependence of the hairpin insertion. 
Indeed, the time dependence of the propagators 
is remarkably well described by the pure double- 
pole formula (4) over the entire observable range 
of time separations. An example for the 16 3 x 32 
lattice at n = .1680 is shown in Fig. 2. 

By fitting the propagators to (|J) , and dividing 
out apppropriate factors obtained from the va- 
lence pion propagator, we obtain a value for mo at 
each k. The 16 3 results are consistent with those 
reported last year Q . The finite volume increase 
observed in the 12 3 data [|| is found to be largely 
an effect of having more nearby visible poles on 
the smaller box. After the MQA shift, the results 
on the two box sizes are within a standard devi- 
ation of each other. The MQA analysis of / w is 
not yet completed, but using previous results for 



/tt, calculated on the 12 3 x 24 configurations, we 
obtain an effective chiral log parameter of 5 ~ .04 
at k = .168. The full results and comparison with 
other work will be presented elsewhere. 

As a byproduct of the hairpin propagator cal- 
culation, we may calculate the integrated pseu- 
doscalar charge density Q§ = j ip^ip d 4 x on 
each lattice in the ensemble. As first suggested 
by Smit and Vink ||, this provides a fermionic 
method for determining the winding number v 
of a gauge configuration, using the integrated 
anomaly equation, which gives v = lim„ w o fnQh- 
Using the allsource quark propagators, improved 
by the MQA pole shifting procedure, we have 
employed the Smit- Vink method to calculate the 
topological susceptibility Xt — (v 2 )/V (where V 
is the lattice 4- volume) for both the 12 3 x 24 
and 16 3 x 32 ensembles. The results show only 
a mild k dependence and can be sensibly ex- 
trapolated to zero quark mass. Using the scale 
a -1 = 1.15 GeV taken from charmonium, we ob- 
tain X t = 12.8 ± 1.4 x 10~ 4 GeV 4 for the 12 3 x 24 
lattice, and Xt = 10.8 ± 1.2 x 10~ 4 GeV 4 for 
the 16 3 x 32 lattice. This can be compared with 
11.5 x 10~ 4 GeV 4 from the WV formula (1). 
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